All multipartite Bell correlation inequalities for two dichotomic observables per site 



R. F. Werner* and M. M. Wolff 

Institut fur Mathematische Physik, TU Braunschweig, 
Mendelssohnstr.3, 38106 Braunschweig, Germany. 
(February 1, 2008) 



O 
O 

(N 
X> 

in 



> 

(N 
O 
(N 
O 



^ • 

9 L,: 
^— > ■ 

c ■ 

=3 : 
cr 



X 



We construct a set of 2 2 independent Bell correlation in- 
equalities for n-partite systems with two dichotomic observ- 
ables each, which is complete in the sense that the inequalities 
are satisfied if and only if the correlations considered allow a 
local classical model. All these inequalities can be summa- 
rized in a single, albeit non-linear inequality. We show that 
quantum correlations satisfy this condition provided the state 
has positive partial transpose with respect to any grouping of 
the n systems into two subsystems. We also provide an effi- 
cient algorithm for finding the maximal quantum mechanical 
violation of each inequality, and show that the maximum is 
always attained for the generalized GHZ state. 

03.65.Bz, 03.67.-a 



I. INTRODUCTION 

Entanglement has not only been a key issue in the on- 
going debate about the foundations of quantum mechan- 
ics, started by Einstein, Podolsky and Rosen in 1935 |Q. 
It also plays a crucial role in the young field of quantum 
information theory. Here entangled states are one of the 
basic ingredients of quantum information processing, due 
to their role as a resource in quantum key distribution, 
super dense coding, quantum teleportation and quan- 
tum error correction (cf. Q). Although general struc- 
tural knowledge about entanglement has improved dra- 
matically in the last few years, there are still many open 
problems. For example, there is still no efficient general 
method to decide whether a given state is entangled or 
not. 

The first, and for a long time also the only mathe- 
matically sharp criteria for entanglement were the Bell 
inequalities They provided the first possibility to 
distinguish experimentally between quantum mechanical 
predictions and those of local realistic models. But al- 
though Bell inequalities have been known for more than 
thirty years [Q, our knowledge about the precise bor- 
der between the classical and quantum mechanical ac- 
cessible region is still mainly restricted to the simplest 
non-trivial cases. Best known is the case of two sites, 
at each of which two dichotomic observables are chosen. 
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This is characterized completely by the Clauser-Hornc- 
Shimony-Holt (CHSH) version of Bell's inequalities [0, 
in the sense that the inequalities are satisfied if and only 
if a local classical model exists f|. Finding a complete 
set of linear inequalities in more complicated situations 
(more sites, more observables, more outcomes) turns out 
to be a very difficult problem in the sense of computa- 
tional complexity M. There is only very little knowl- 
edge about Bell type inequalities beyond the CHSH case 
|p[-|l2|. Though numerical studies yield a large number 
of inequalities p3| , for most of them it is neither known 
by how much they can be violated in quantum theory 
nor is there a general characterization admitting further 
investigations. 

We were therefore quite surprised ourselves at finding 
an infinite sequence of multipartite correlation settings 
for which we could develop the theory to be as explicit 
and complete as in the CHSH case. Our setting general- 
izes the CHSH-setting to an arbitrary number n rather 
than two different sites, but retains the constraints of just 
two observables per site with just two outcomes each. 
Thus each of the n-participants has the choice of two ob- 
servables, each of which can take the values +1 or — 1. 
For any choice of observables we then consider the ex- 
pectation value of the product of all n signs (a "full" 
correlation function). A Bell inequality is a linear con- 
straint on the set of all such expectations, which is valid 
whenever the correlations can be obtained from a local 
classical model, and which cannot be written as a con- 
vex combination of other such constraints. Examples are 
the CHSH inequality || for n — 2 and their generaliza- 
tions going back to Mermin and others f|-|ll| leading to 
a single inequality for arbitrary n. 

We remark that this problem setting could be general- 
ized to include the expectations not only of the product 
of all n signs, but also the products of subsets of signs ( 
"restricted" correlation functions). These data would be 
sufficient to reconstruct the full joint probability distri- 
butions of signs for all choices of observables. However, 
most of the derivations in this paper do not generalize 
to this setting, and it is not yet clear which statements 
would still be valid (maybe with a different proof). When 
we talk of the existence of a classical model, however, it 
is understood that such a model would also determine 
all restricted correlation functions. The omission of re- 
stricted correlation functions from our setting only means 
that we do not consider constraints depending on them. 

For this class of multipartite correlations we obtained 
the following results: 
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• We construct a set of 2 2 ™ Bell inequalities, and 
show its completeness: the correlations considered 
allow a local classical model if and only if all these 
inequalities are satisfied (Sec. |lj]). 

• The convex set of collections of classical correlation 
functions is a 2™-dimensional hyper-octahedron, 
which can be described alternatively by a single 
nonlinear inequality (Sec. III). 



• We discuss the symmetries connecting different 
inequalities and develop a construction scheme, 
which yields all 2 2 equalities by successive sub- 
stitutions into the CHSH inequality (Sec. [Tv| ). 

• We reduce the computation of the maximal quan- 
tum violations of each Bell inequality to a simple 
variational problem with just one free variable per 
site. The maxima are already attained in qubit 
systems, more specifically for the n-party general- 
ization of the GHZ state |Q, with a choice of ob- 
servables depending on the inequality under con- 
sideration (SecM). 

• We extend this method to a characterization of the 
convex body of quantum mechanically attainable 
correlation functions in terms of its extreme points, 
which are also found in the generalized GHZ state. 
These results are analogous to those of Tsirelson 
[p^| , ^6| for the bipartite case. 

• We characterize the Mermin inequality as that Bell 
inequality, which can be violated by the widest 
margin in quantum theory. 



• Sec. VI settles the relationship between the corre- 
lation Bell inequalities and another important en- 
tanglement property. We show that for states hav- 
ing positive partial transposes with respect to all 
their subsystems, all 2 2 inequalities are satisfied, 
so the correlations in such quantum states can be 
explained in the context of a local realistic model. 
This extends our earlier result Jl7t for Mermin's in- 
equalities, and is further supporting evidence for a 
recent conjecture by Peres [0, namely that posi- 
tivity of partial transposes should generally imply 
the existence of local realistic models. 

In the appendix we will discuss some of the general 
results obtained in the sections III , |v|, [v| in more detail 
for the special cases n = 3, 4. 



II. BELL'S INEQUALITIES AND CONVEX 
GEOMETRY 

Before entering the discussion of Bell inequalities in our 
special context it is useful to recall some geometric struc- 
tures of the general problem and basic facts concerning 
the duality of convex polytopes. 



Consider a system decomposed into n independent 
subsystems. Suppose further that on each of these sub- 
systems one out of m v- valued observables is measured. 
Thus each of the m n different experimental setups may 
lead to v n different outcomes, so that the raw experimen- 
tal data are made up of {mv) n probabilities. These num- 
bers form a vector £ lying in a space of dimension (mv) n 
(minus a few for normalization constraints). Classically, 
in a local realistic model, £ would be generated by spec- 
ifying probabilities for each classical configuration, i.e., 
for every assignment of one of the v values to each of the 
nm observables. Here the "local" character of the the- 
ory is expressed by the property that the assignment of 
a value to an observable at site k does not depend on the 
observables chosen at other sites. Every configuration 
c also represents a possible classical (ideally prepared) 
state, and hence a vector e c of probabilities. The classi- 
cal accessible region, which we will denote by Q, is thus 
the convex hull of i/"™ 1 ) explicitly known extreme points. 
Even though the number of configurations is large, it is 
finite, hence f2 is a polytope. 

Like every compact convex set, f2 is the intersection of 
all half spaces containing it. A half space is completely 
characterized by a linear inequality, so we must look for 
vectors (3 such that ((3, £) < 1 for all £ G fl. Since this 
property can be checked on the extreme points e c we must 
look at the convex set 



B={f3\Vc:((3,e c )<l}, 



(1) 



also known as the polar of {e c }. For each j3 £ B the 
inequality (/?, £) < 1 is thus a necessary condition for 
£ € £1. Moreover, the Bipolar Theorem []19| says that the 
collection of all these inequalities is also sufficient. 

Luckily, the inequalities are not all independent, since 
the inequality for a convex combination (3 = ^Xifii, with 
(3i G B already follows from the inequalities for the Pi. 
It therefore suffices to take only the extreme points of 
B. For a polytope this has a very intuitive geometrical 
interpretation: the half spaces determined by extreme 
points touch Q in a face of maximal dimension. Moreover, 
there are only finitely many such maximal faces, which 
is to say that B is also a polytope. 

The task of finding all Bell inequalities is therefore a 
special instance of a standard problem in convex geome- 
try, known as the hull problem: given the extreme points 
{e c } of a polytope fi, find its maximal faces or, equiva- 
lently, the extreme points of its polar. 

The duality between B and O is a generalization of 
the duality between regular platonic solids, under which 
dodecahedron and the icosahedron, as well as the octa- 
hedron and the cube are polars of each other. A gen- 
eralized (d-dimensional) octahedron is the unit sphere 
in a sequence space ^ 1 ({l,...,d}). Its polar is the unit 
sphere in the dual Banach space £°°({1, . . . ,d}), i.e., a 
rf-dimensional hypercube. This is precisely the situation 
we will find for the classically accessible region considered 
in this paper, where d = 2". 
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The first to consider the construction of a complete set 
of Bell type inequalities as a problem in convex geometry 
apparently was M. Froissart pofl . Unfortunately how- 
ever, a general solution for all (n, m, v) is highly unlikely 
to exist. To find some extreme points of ([l]) is not so diffi- 
cult, but algorithms providing the complete set are likely 
to run into serious growth problems already for very small 
(n,m,v). In fact, there is a theorem by Pitowsky jjj to 
the effect that, in a closely related problem, finding all 
inequalities would also solve some known hard problems 
in computational complexity (this is in fact strongly con- 
nected with the notorious NP = P resp. NP = coNP 
questions). Pitowsky and Svozil |l3fl have recently per- 
formed an extensive numerical search for n = 3, and 
published their result, the coefficients of 53856 inequal- 
ities on their website. Unfortunately, there is not much 
generalizable insight coming out of this kind of work, but 
it is nice to see what can be done in this hard numeri- 
cal problem. For further problems and partial results in 
this genre we refer to the problem page |T^ | on our own 
website. 

In what follows we will restrict to the case (n, m, v) = 
(n, 2, 2) and "full" correlation functions in the sense de- 
scribed in the introduction. 



III. ALL BELL CORRELATION INEQUALITIES 
A. Basic notation 

Talking about Bell inequalities one usually has in mind 
inequalities of the Clauser-Horne-Shimony-Holt form 0. 
These inequalities refer to correlation experiments, in 
which each of two parties has the choice of two ±1 valued 
observables to be measured, i.e., (n, m, v) = (2, 2, 2). Fo- 
cusing only on full correlation functions for multi-particle 
generalizations of such systems ((n,m,v) = (n, 2,2), n 
fixed arbitrarily) the raw experimental data are 2" ex- 
pectation values, each corresponding to a different exper- 
imental setup. Each setup is labeled by the choice of ob- 
servables at each site. We parameterize these choices by 
binary variables s k G {0, 1} so that Sk indicates the choice 
of the ±l-valued observable Ak(sk) at site k. Each full 
correlation function is thus the expectation of a product 
Y\ k Ak{sk), and is labeled by a bit string s — (si, . . . , s„). 

We will consider these expectations as the components 
£ (s) of a vector £ in a 2™-dimensional space. Then any 
Bell inequality is of the form 

5>(«)*00<1, ( 2 ) 

s 

where we have normalized the coefficients (3 so that the 
maximal classical value is 1, in accordance with the def- 
inition of polars in Sec. [n]. The linear combination in 
Eq. (||) can also be computed under the expectation 
value, so that this inequality can be stated as an upper 
bound on the expectation of 



B = Y,P( s )I[Msk). (3) 

s fc=l 

We call such expressions Bell polynomials. They can be 
used directly in the quantum case, where all variables 
Ak(sk) are substituted by operators with —1 < Ak(sk) < 
1, acting in the Hilbert space of the fc-th site, and the 
product is taken as the tensor product. It is often useful 
to consider these polynomials rather than the set of coef- 
ficients, because often many coefficients are zero, and we 
can sometimes simplify a polynomial algebraically (e.g., 
by factorization), even though this may not be apparent 
from the coefficients. 

Two convex sets in the real 2"-dimensional vector 
space are the subject of our investigation: firstly, the 
polytope fl of correlation vectors £ coming from local 
classical models, and secondly the set Q D 57 of such 
vectors arising from quantum models. Q will be charac- 
terized in terms of Bell inequalities in this section, Q will 
be considered in Sec. M. 



B. Construction and completeness 

In a local classical model every observable A k (sk) is 
a random variable in its own right, i.e., it is a func- 
tion of the "hidden variable" which does not depend on 
the choices si of observables at other sites i ^ k. A 
model must assign probabilities to any collection of val- 
ues for these observables, i.e., to each classical configura- 
tion. Since the extremal choices of such probabilities just 
assign probability 1 to one configuration and zero prob- 
ability to all others, the extreme points of f2 are simply 
labeled by the configurations. 

One configuration c is the choice of Cfc(sfc) G {— 1, 1} 
for all k and Sfc. Clearly, there are 2 2 ™ such configura- 
tions. The corresponding correlation vector £ = e c has 
components 

n 

tc{s) = ~[[c k {s k ) . (4) 

fc=i 

Since we only consider full correlation functions (and not 
restricted ones, see the introduction), different classical 
configurations may give the same extreme point e c . For 
example, we may choose two different sites, and change 
the values of all Ck(sk) at these sites simultaneously. 
Then in Eq. (Q) the sign changes cancel for all s. This is 
also apparent from the factorization 

(n \ n 

k=l I 1=1 

in which the first factor is just an s-independent sign, and 
in the second factor it suffices to choose configurations 
with Cfc(0) = 1. Thus we can write Ck(sk) — {— l) Sfcrfc 
with r k e {0,1}. Then 
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e c (s) = ±(-l) 



E 

k=l 



rkSk, 



(0) 



where the extreme points are now labeled uniquely by 
the bit string r = {r\, . . . ,r n ) and the overall sign. This 
leaves us with exactly 2 n+1 extreme points of f2. 

Our task is now to find the extremal linear inequalities 
/?, characterizing this set, i.e., the extreme points of B 
from Eq. (|l]). The bipartite case was indeed completely 
analyzed by Fine Q , who showed that there are only two 
classes of inequalities: one is trivial in the sense, that it 
just requires correlations to be in [— 1,+1], and the the 
second consists of the CHSH type inequalities, for which 
the prototype is (3 — (ftfl' - 2)' ^ construction of 
some Bell type inequalities for arbitrary n was first pro- 
posed by Mermin || and further developed by Ardehali 
[§, Klyshko @ and Gisin 0. 

We will now find all extremal solutions (3 to the set of 
inequalities 



l<X>(*)(-l)< r '*> <1, 



(7) 



where r G {0,1}™ runs over all bit strings characterizing 
the configurations. Suppose that of these 2™ inequali- 
ties p < 2 n fixed ones are "tight" in the sense that the 
sum takes one of the extreme values ±1. This will be 
consistent with a plane (affine manifold) of vectors (3 of 
dimension at least 2 n — p. We can now construct convex 
decompositions of (3 in an open neighborhood of (3 in this 
plane, since each one of remaining sums is continuous in 
(3, and there is a finite margin before another inequal- 
ity becomes violated. This contradicts extremality, so 
we conclude that the inequality must be tight for all r. 
Thus we have 2" signs f(r) e {+1,-1} with 



(r,s) 



f(r), 



(8) 



For examples of this numbering, see the appendix. The 
converse procedure is similar. For example, the Math- 
ematica package available from our website |12| finds 
that Mermin's inequality for n — 6 has the number 
1 692 930 046 964 590 721. 



C. Structure of the classical region 

From the previous section it is clear that the classi- 
cal region O is a polytope in d = 2™ dimensions with 
2d extreme points and 2 d maximal faces. This suggests 
that fl should be a hyper-octahedron, whose polar B is a 
hyper-cube. Indeed from the parametrization of the in- 
equalities by d values f(r) = ±1, the latter statement is 
rather obvious. That fl is an octahedron is not so appar- 
ent in the coordinates labeled by s as above. However, 
we can choose a basis transformation making this geo- 
metric identification of f2 more obvious. The necessary 
transformation is, of course, just the Fourier transform. 
With the notation 



|(r) = 2-"£(-l)<^(s) 



(10) 



we can summarize the findings of the previous section by 
saying that £ £ fl if and only if 



V/G{-1,1} 2 ": £/(r)|(r)<l. 



(11) 



The expression in (|llj) reaches its maximum with respect 
to /, if f(r) is just the sign of £(r). Therefore, the whole 
set of 2 2 linear inequalities (or the statement £ S 51) is 
equivalent to the single non-linear inequality 



Ek»i 



< 1. 



(12) 



Now we can read Equation (||) as a Fourier transform 
with respect to the group of n-tuples of {0, 1} with addi- 
tion modulo 2. Therefore, we easily obtain the entire set 
of extremal (3 by applying the inverse transformation to 
the set of vectors / € { — 1, l} 2 : 



(r,s) 



(9) 



These are the coefficients of the complete set of 2 2 " ex- 
tremal Bell inequalities specifying the range of expecta- 
tions of full correlation functions for any local realistic 
model. 

The inequalities constructed in this way have a natu- 
ral numbering, defined by the following procedure: For 
any number between to 2 2 — 1, write the binary ex- 
pansion with "digits" ±1 to get /, and perform the in- 
verse Fourier transform (^|). From (3 compute the poly- 



nomial (t 
equality, 



), which is often the best form of writing the in- 
Decause one can apply algebraic simplifications. 



Obviously, this nonlinear inequality is nothing but the 
characterization of the hyper-octahedron in 2™ dimen- 
sions as the unit sphere of the Banach space t 1 . 

From this simple characterization of fl it might seem 
that our problem is essentially trivial. However, the vast 
symmetry group of £1, which includes among other trans- 
formations the set of (2")! permutations of the coordi- 
nates is misleading, because these are not really symme- 
tries of the underlying problem of finding all correlations 
within a classical model. This is apparent from the obser- 
vation that the Bell polynomials associated with the ex- 
treme points may look quite different algebraically. That 
is, the 2™ dimensions are not really equivalent, but carry 
some structure coming from the division of the system 
into n sites. This is even more obvious when looking at 
the set of quantum correlations, which has a much lower 
symmetry. 

Nevertheless, the underlying problem has a large sym- 
metry group, which will be studied in the next section. 
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IV. SYMMETRIES AND SUBSTITUTIONS 

Browsing through the complete set of linear correla- 
tion inequalities one quickly gets the feeling that there 
are many rather similar ones, and also some inequali- 
ties which can be obtained in a rather trivial way (e.g., 
as a product) from lower order ones. In this section we 
will describe the grouping of the inequalities into "essen- 
tially different ones" , and also how they can be obtained 
by an efficient construction for composing higher order 
inequalities from lower order ones. Both ways of struc- 
turing the set of inequalities make sense for more general 
cases (n,m,v) (see Section ||), but for the moment we 
only apply them to our restricted class. 

A. Symmetry Group 

Some symmetries acting on Bell inequalities are ob- 
vious and, in fact, present in any problem of this type, 
involving any number of outcomes and observables. The 
basic symmetries leading to equivalent inequalities are: 

(i) . Changing the labeling of the observables at each 

site. 

(ii) . Changing the names of the outcomes of each ob- 

servable. 

(iii) . Permuting subsystems. 

Since we have two observables per site, there are 2™ ways 
of swapping the labels of observables at each site. Swap- 
ping the ±1 outcomes of an observable Ak(sk) at site 
k results in a sign in all correlation functions involving 
this observable. We have already utilized the fact that 
swapping both Ak(0) and Afe(l) only results in an overall 
sign, so it is enough to consider sign changes for Ak(l) 
only. Clearly, there are 2™ such sign changes. Expressed 
in terms of the function / these transformations amount 
to 

/(r )^(_l)<*o,r> /(r + ro)5 (13) 

where r, so ; r o all lie m {0,1}™, and ro and sq are the 
parameters describing the sign changes and observable 
swaps, respectively. Together with the global sign change 
and the n! permutations we thus find the group G of 
symmetry transformations in our case to have the order 

\G\ = n\ 2 2n+1 . (14) 

The orbit of a given inequality is defined as the set of 
all the inequalities generated from it by symmetry trans- 
formations. The number of elements in an orbit is \G\, 
divided by the order of the group of symmetries leav- 
ing an element of the orbit invariant. The number of 
different orbits is the number of "essentially different" 
inequalities. Obviously, (|lj) is an upper bound on the 



number of elements in each orbit. Since the union of all 
orbits is the set of all inequalities, this leads to a lower 
bound on the number of essentially different inequalities. 

Note that \G\ increases much more slowly than 2 2 , 
the total number of extremal inequalities. Therefore, for 
large n the classification up to symmetry hardly reduces 
the number of cases. Explicitly, we find: 



n 


inequalities 


\G\ 


orbits 


2 


16 


64 


2 


3 


256 


768 


5 


4 


65 536 


12 288 


39 


5 


4 294 967 296 


245 760 


> 17 476 



For n up to 4, the number of orbits was obtained explic- 
itly. However, for n > 5 the lower bound on the number 
of orbits makes it clear that listing all essentially differ- 
ent inequalities is not going to be useful. More detailed 
results up to n = 4 will be shown in the appendix. 

B. Generating new inequalities by substitution 

A simple way of generating inequalities for higher n 
is to partition the n sites into two subsets of sizes n\ 
and n 2 = n — Ui and to take arbitrary Bell polynomials 
for rii and n 2 sites, appropriately rename the variables, 
and to multiply the two expressions. For example, the 
polynomial 

i(aibi + a x b 2 + a 2 bi - b 1 b 2 )ci (15) 

is obtained by multiplying a CHSH polynomial for the 
first two sites with the trivial polynomial "ci" on the 
third (note that for the sake of clarity we have substituted 
A 1 (0),A 2 (1) with ai,b 2 etc.). It is clear that this gives 
an extremal Bell inequality for three sites. 

This procedure can be generalized considerably by not- 
ing that the product operation corresponds to the trivial 
two site Bell polynomial "ai&i", but nothing restricts us 
to using a trivial expression here. So in general, con- 
sider a partition of the sites into K subsets of sizes rik, 
X^jtLi n k = n - Then pick an extremal Bell polynomial for 
K sites, written out in variables Ai(0), Ai(l), . . . , Ak(X)- 
Now substitute for each Ak(sk) an extremal Bell polyno- 
mial for rik sites. We claim that the resulting polynomial 
in n variables is an extremal Bell polynomial. 

Indeed, if we substitute for each of the variables either 
+1 or —1, we will get Ak(sk) = ±1 for each k, Sfc, because 
we substituted extremal Bell polynomials. But then the 
same argument on the level of K sites shows that the 
value will be ±1. 

We will say that a Bell polynomial is elementary, if 
it cannot be obtained by substitution from lower order 
polynomials. Obviously, if an inequality is elementary, 
so is its entire orbit. Clearly, the CHSH inequality is ele- 
mentary. Moreover, it is known that it is a good tool for 
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generating higher order inequalities by substitution: one 
of the constructions jl(]Jll]] of the Mermin's inequalities 
is based on this idea. But in view of the rapid increase of 
the double exponential one might think that there must 
be many more elementary inequalities. However, we have 
the following result: 

Proposition. The CHSH-inequality is the only ele- 
mentary Bell inequality in the class we consider, i.e., all 
these inequalities for n > 2 can be constructed by succes- 
sive substitutions into the CHSH-inequality. 

It is an interesting open problem, whether this state- 
ment holds for other families of Bell inequalities, e.g., the 
one tabulated in | fL3|| . 

We start the proof on the level of vectors / € {— 1, l} 2 
parameterizing an arbitrary extremal Bell inequality for 
n sites. We decompose the system into a partition of 
K = 2 subsets of size n — 1 , 1 and rewrite 



f(n, r n _i, r n ) = f(r, 0)S rnt0 + f(r, l)5 rn 



(16) 



The respective coefficients /3(s) of the n-site inequality 
are then obtained via Fourier transformation according 
to Eq. (g), and we get 

j3{s) = 2-"]T/(r)(-l)^> 

r 



/3 (5} + (-l) s "/9i(5) 



(17) 



where [3)~ (s) are coefficients for extremal Bell inequalities 
for n — 1 sites. If we now add the respective observables 
Ak(sk) and write out the corresponding Bell polynomial 



s=$>(s)n^fc(s fe ) 

fc=l 

- -S [a„(0) + A n (l) 



-Bi 



A n {0)-A n (l) , (18) 



we immediately see, that this is just a CHSH polynomial, 
where the observables of one site have been substituted 
by Bell polynomials Bq and B\ for n — 1 sites. 



V. QUANTUM VIOLATIONS 

Provided with a huge number of Bell type inequalities 
we now go beyond the classical accessible region. The 
first question to arise is of course whether or not and to 
what extent quantum systems can violate these inequal- 
ities. To answer this question we will first provide an 
effective variational method for computing the maximal 
quantum violations and show, that they are bounded by 
those obtained for Mermin's inequalities. In the following 



two subsections we will then briefly discuss the structure 
of the underlying quantum domain, and prove that the 
generalized GHZ state maximally violates any of the cor- 
relation inequalities. 



A. Obtaining the maximal violations 

In order to compute the maximal quantum violation 
of any correlation inequality we have to vary over one 
density operator p on a tensor product of n factors, and 
two operators in each factor. Assuming all tensor factors 
to have dimension d, this means d 2n parameters for the 
density operator and 2nd 2 for the observables. Hence the 
numerical solution of this variational problem is not fea- 
sible, except for the most trivial cases (and even impos- 
sible, because d is, in principle, a free parameter). Fortu- 
nately, however, it turns out that computing the overall 
maximum is much easier than computing the maximal vi- 
olation for a fixed state: we will reduce the computation 
to a variational formula in just n variables. 

First we have to recall some basic notions. In quantum 
mechanics expectations of ±l-valued observables are de- 
scribed by Hermitian operators Ak(-Sk) with spectrum in 
[—1, +1]. Since we are only interested in maximal corre- 
lations, we may as well take the observables extremal in 
the convex set of Hermitian operators with — 1 < A < 1, 
i.e., we may assume the observables to be unitary and 
thus A 2 = 1. 

The general form of a Bell inequality for an n-partite 
quantum system, which is characterized by a density op- 
erator p, is then 



n 

tr(pB): = tr[p^/3( S )(g) 



A k {s k ) 



fe=i 



< 1, 



(19) 



where we will refer to B as the Bell operator, which is just 
the quantum counterpart of the Bell polynomial defined 
in Eq.(H). Of course, every expectation value ( |l9| ) larger 
than 1 is called a violation of Bell's inequality. 

In order to derive the maximal quantum violation, 
which is nothing but the operator norm of the Bell oper- 
ator, we first define another operator C by 



C:=B(g)A k (0) 



)C s k k 



(20) 



fc=i 



fc=i 



where we have set C% = Ak(l)Ak(0), and C° = I. Since 
the Ck are commuting unitary operators, all summands 
of C can be diagonalizcd simultaneously, and the eigen- 
vectors of C are tensor products of eigenvectors of the 
Cfe . Every eigenvalue 7 of C is therefore of the form 



s fc=l 



(21) 



where jk is an eigenvalue of Ck ■ It is clear from the above 
remarks that C commutes with its adjoint, so ||C|| is just 



G 



the modulus of the largest eigenvalue. Now we utilize the 
fact that \\B*B\\ = \\C*C\\, i.e., ||B|| = ||C|| and obtain 



\B\\ 



sup 



e^ik* 

s k=l 



(22) 



where each 7^ runs over the eigenvalues of Ck = 
Ak(l)Ak(0). This formula allows us to compute the 
largest expectation tr(pB) for fixed real coefficients (3 
(coming from a Bell inequality or not) and a fixed choice 
of observables Ak(sk), but with p chosen without further 
constraints to maximize the expectation. 

What we are now interested in is the maximum also 
with respect to the A&(sfc). Since formula J22] ) depends 
only on the eigenvalues 7& this will be given by the same 
expression, but with 7^ running not just over the eigen- 
values of a particular operator C k but over all 7^ which 
can be eigenvalues of products of unitary and hermitian 
operators. Since such a product is again unitary, we have 
1 7k I = 1. Moreover, as is easily seen in 2 x 2 examples, 
this is the onl y con straint in any Hilbert space dimension 
(see also Sec. VD). Hence for any choice of real coeffi- 
cients f3(s) and observables — 1 < Afc(sfc) < I, we have, 
as the best possible bound: 



n 

tr[p]T/3( S )(gM fe ( Sfc ) 



k=i 



< sup E^n 

fTfc> a fe=l 



St 

Ik 



(23) 



where the supremum runs over all {7i,---,7n} with 
|7fe| = 1. Moreover, the bound does not change with 
Hilbert space dimension, as long as all factors are non- 
trivial. 

A more detailed discussion of quantum violations uti- 
lizing Eq. (^||) for the special cases n — 3, 4 can be found 
in the appendix. 



B. Mermin's inequalities and the overall maximum 

Asking for the overall maximal quantum violation we 
may additionally vary over the set of inequalities. In 
utilizing the result obtained in the previous section, 
we are able to express the norm of a Bell operator in 
terms of lower order Bell operators. Moreover it suf- 
fices to consider qubit systems and we may therefore set 
Ak(sk) — afe(sfc)(?, where o is the vector of Pauli matrices 
and Sk(sk) is a normalized vector in R 3 . 

Squaring Eq.([l8|) this leads to 



B z 



B 2 B 2 

® [1 + a n (0)a n (l)] + -± ® [1 - a„(0)a n (l)] 



+ [Bo, Si] <g> -[a„(0) x a„(l)]a. 



(24) 



Without loss of generality we now assume that ||i?o|| < 
|| Bill and estimate by induction 



This bound is indeed saturated by the set of inequal- 
ities going back to Mermin [js|-[il"||, which thus provide 
the overall maximal quantum violation. In fact, we will 
show, that the converse is also true, so that we have the 
following 

Proposition. The orbit corresponding to Mermin's in- 
equality is the only one for which the maximal violation 
2 J ^ — is attained. 

Before we continue proving the claimed uniqueness, we 
emphasize, however, that this does in general not imply, 
that for a fixed quantum state Mermin's inequality is 
more strongly violated than any other. 

We begin our proof with noting that the maximal norm 
of the Bell operator in Eq.(p4[) requires orthogonality 
of the observables, such that the respective phases in 
Eq.(^|[) have to be ±i. Without loss of generality we 
can thereby restrict to the case +i since the remaining 
sings just correspond to a transformation between two 
inequalities of the same orbit according to (|l3|). Hence, 
Eq.(||) leads to 



\B n 



2-*|E wo ii^- 1 )* 

r,s k—1 

n 

2- n \j2f max (r)l[[l + i(-ir] 

r k=l 
n 

2-*|E/m«(r) Y[[^ 2rk) ] 

r k=l 



(26) 



If we now want B max to saturate the bound in ( |25| ) , then 
following Eq.(]26]) we are left with four possible choice 

for f max , like f ma x{r) = 1 for (-i)^k Vk = l J and 
fmax{r) = — 1 otherwise. Since these four inequalities 
again belong to the same orbit, the correlation inequal- 
ity leading to the overall maximal quantum violation is 
indeed uniquely determined (up to equivalence transfor- 
mations within one orbit). 



C. Structure of the quantum domain 

In the same manner as we did for the classical case we 
may ask for the structure of the region in the space of 
correlations, which is accessible within the framework of 
quantum mechanics. One of the first to investigate this 
question in more detail apparently was Tsirelson Jl5| , ^6| , 
while studying quantum generalizations of Bell's inequal- 
ities. 

Let us begin with defining the quantum counterpart of 
the classical accessible region f2, introduced in Sec. ||: 



|B|| 2 = ||B 2 || < 2||Bi|| 2 < 2"- 1 . 



(25) 



6 = tr 



pQ<)A k {s k ] 

fc=i 



C IT 



(27) 
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where {A k } are suitable observables and p is a quantum 
state in arbitrary dimension. The structure of Q is much 
more complicated than that of il C Q. In particular, it 
is not a polytope. Nevertheless, we can explicitly param- 
eterize its extreme points. For the sake of completeness 
we will first prove convexity of Q, although this follows 
closely the work of Tsirclson @]. 

Consider a convex combination of vectors in Q 



(28) 



and an associated Hilbert space 

n 

W = (g)H fe =(g)0?4 Q) - 0(g)«i afe) . (29) 



fe=l fe a 



ai...a n fc 



Then with p = (J) Q A^p , which is a density operator 
acting on the "diagonal subspace" 



®(8)Hi a) cH, 



(30) 



a fc 



and A/^Sfc) = Q Aj^(sfc) we are given a state and ob- 
servables such, that the convex combination in ( p8| ) is 
indeed a proper element of Q. Hence, Q is convex. 

Now let us return to the result obtained in the previous 
subsection. Following Eq. 
quantum violation of an ar 



we can write the maximal 
}itrary inequality (3 as 



sup Y,^ s M^U^ k ) 

a V fe=l 

= sup ^2/3(s)£ s (<p ,---,<Pn), 



(31) 

(32) 



where we have set £ s (<po, . . ■ , <£„) = cos(<p + Sfc ^fc^fe)- 
Now by the Bipolar Theorem |19] the convex set Q is just 
given by the convex hull of these vectors: 

Q = co{£(ip ,...,ip n )}. (33) 



D. Generalized GHZ states 

It is a well known fact, that the generalized GHZ state 
defined by 



1 



\^ GHZ ) = _(|00...0) + |11...1» 



(34) 



maximally violates Mermin's inequalities pffl . Astonish- 
ingly this is also true for any other of the 2r correlation 
inequalities: 

Proposition. Any extreme point of the convex set of 
quantum correlation functions as defined in Eq. $?\) is 



already obtained for the generalized GHZ state. In par- 
ticular, this implies that GHZ states maximally violate 
any of the presented correlation inequalities. 

We have to show, that for any set of angles 
{ipo, ■ ■ ■ ,<fn} there are suitable observables such, that 

n 

(^ghzI® A k (s k )\f GHZ ) = cos^o + VfeSfe)- (35) 
fe=i fc 

Therefore we choose observables A k {s k ) = a k (s k )a with 
flfc(O) = (cos a, sin a, 0) 

a fe (l) = (cos((^ fc + a), sin(<p fc + a), 0). (36) 

These observables simply swap the basis vectors provid- 
ing them with an additional phase factor, i.e., 

a k ( s k)# |i) = exp [i{-l)\a + <p k s k )] \j © 1), (37) 

where j = 0, 1 and © means addition modulo 2. Hence, 
for the left hand side of Eq.(|35|) two terms occur, which 
are just complex conjugates of each other, and we get 

n 

(*GHz\®Msk)\*GBz) =5fe{e M V£^ fcSfc }, (38) 
fc=i 

so that it just remains to set a = fo/n. 



VI. STATES WITH POSITIVE PARTIAL 
TRANSPOSES 

The violation of one of the inequalities, which can be 
derived from Eq.Q, is a rather physical entanglement 
criterion, since we can at least in principal decide it ex- 
perimentally by measuring the respective correlations. 
However, the difficulty in doing so is the choice of the 
observables, and optimizing them for a fixed state leads 
in general to a very high dimensional variational problem. 
An entanglement criterion, which is in contrast easy to 
compute, is the partial transpose proposed by Peres in 

13 ■ 

Before we settle the relationship between these two en- 
tanglement criteria, we will briefly recall some basic no- 
tions. 

The partial transpose of an operator on a twofold tensor 
product of Hilbert spaces Hi ® TL2 is defined by 



D 



Ti 



D 



3' 



(39) 



where Cj on the right hand side is the ordinary trans- 
position of matrices with respect to a fixed basis. The 
generalization of this definition to an n-fold tensor prod- 
uct is straight forward, and we will denote the transpo- 
sition of all sites belonging to a set r C {1, . . . , n} by the 
superscript T T . 



8 



Recall further, that a state is called separable or clas- 
sically correlated, if it can be written as a convex com- 
bination of tensor product states - otherwise it is called 
entangled. A necessary condition for separability, which 
also turned out to be sufficient in the case of two qubits 
p2f , but not in general (cf. j23| ) , is the positivity of all 
partial transposes with respect to all subsystems. More- 
over, there is a conjecture by Peres jL8|, that this might 
even imply the existence of a local realistic model. In 
jl7j we showed, that the set of inequalities going back 
to Mermin is indeed fulfilled for states satisfying this 
"ppt" -condition. In the following we will show that this 
implication is not due to a special property of these in- 
equalities, but holds for any Bell type inequality in (19), 
as long as we consider expectations of full n-site correla- 
tions. This leads to the main result of this section: 

Proposition Consider an n-partite quantum system, 
where each of the parties has the choice of two dichotomic 
observables to be measured. Assume further, that the par- 
tial transposes with respect to all subsystems of the corre- 
sponding density operator are again positive semi-definite 
operators. Then all the 2™ correlations can be described 
in the context of a local realistic model. 

In particular, this implies, that if a state is biseparable 
with respect to all partitions, all the inequalities are sat- 
isfied even if there exists no convex decomposition into 
n-fold product states. 

In order to prove this proposition and to derive an 
upper bound for the expectation of the Bell operator, we 
first apply the variance inequality to p Tt and B Tt : 

(tvpB) 2 = (trp^B^f < tv[p T -{B T -) 2 ] 

<tr{p[(B^}. (40) 

Since we suppose that p Tt > Vr this holds for any 
partial transposition, and we may take the average with 
respect to all subsets r, and have therefore to estimate 
the expectation of the operator 



^ E E <8> Msk)A k (s' k ) (g) A k (s' k )A k (s k ) 

T s,s' k£r k$:T 

n , 

= J2P( s M s ')®o{Msk),A k (s' k )} + , (41) 



k=l 



where {■, -} + denotes the anti-commutator. Note that in 
the first line of Eq.(f4l"|) we have rearranged the tensor 
product and made use of 

[A k {s' k ) T A k {s k ) T ] T = A k {s k )A k {s' k ). (42) 

Since A 2 — 1 and s k ,s' k € {0,1} only two different op- 
erators can arise in every tensor factor in ([H|): either 
\{A k {Q) , A k {\)} + or the identity operator. These two 
obviously commute, and we can therefore simultaneously 
diagonalize all the summands. What remains to do, is to 



substantiate our intuition that "if everything commutes, 
then we are in the classical regime" . For this purpose 
note, that eigenvalues of the operator (Elh are of the form 



E/w^nn* 



fe=l 



. s k 



}■ 



(43) 



for suitable — 1 < Xk < 1- But since we can always find 
classical observables C with correlations 



(C fc (0)C fc (l)) = X k, 



(44) 



we are able to construct a system, which is classical in the 
sense that it may be described in the context of classical 
probability theory, such, that (^) is the expectation of 
the square of the respective Bell polynomial. However, 
due to the defining properties of the Bell inequalities, 
this is indeed bounded by unity, which proves our claim, 
that all the considered Bell inequalities are satisfied for 
states having positive partial transposes with respect to 
all their subsystems. 



VII. CONCLUSION 

We provided two approaches for constructing the en- 
tire set of multipartite correlation Bell inequalities for 
two dichotomic observables per site: the Fourier trans- 
formation of a 2™-digit binary number and nesting CHSH 
inequalities. This set of inequalities led us to a single non- 
linear inequality, which detects the existence of a local 
classical model with respect to the considered correla- 
tions. We were able to simplify the variational problem 
of obtaining the maximal quantum violation of the linear 
correlation inequalities, in particular showing, that these 
arc attained for generalized GHZ states, and proved, that 
"ppt states" satisfy all these 2 2 " inequalities. 

One crucial assumptions was, that each site has the 
only choice of two dichotomic observables to be mea- 
sured. Permitting more observables per site, more out- 
comes per observable or even the choice of " not mea- 
suring", i.e., including restricted correlation functions, 
would lead to non-commuting terms, and most of the ar- 
guments would fail. So this is obviously a starting point 
for further investigations. In particular, one may think of 
applying the mechanism of substitution (Sec.|V|) in order 
to derive new classes of Bell inequalities. 

Another open question concerns the hierarchy of the 
inequalities with respect to their quantum violations. 
That is, if a given inequality is violated for a fixed quan- 
tum state, is there a set of inequivalent inequalities, 
which have to be violated as well? 

Finally, we want to mention, that there is recent work 
by Scarani and Gisin j24j] pointing out, that there might 
be a close relation between the quantum violation of mul- 
tipartite Bell inequalities and the security of n-partner 
quantum communication. 
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APPENDIX: 

Recently more and more attention has turned to tri- 
and four-partite states, especially to symmetric states as 
laboratories for multipartite entanglement (cf. [f25|~p7f). 
Therefore we will provide the complete set of Bell in- 
equalities for these cases in a more explicit form and ad- 
ditionally give the maximal quantum violations, which 
we have numerically [2S[] obtained utilizing the method 
presented in Sec]y|. 



1. Inequalities for three sites 

For n — 3 Eq.(|^) leads to the five essentially differ- 
ent Bell polynomials (for the sake of legibility we again 
substitute Ai(0), ^2(1) with 01,62 etc.): 



2 Yl °* 6 - 



ic: n 



k,L: 



2 M & i 



Oi&iCi 

02(^1 - 62)] Cl 



- [ai6i(ci + c 2 ) - a 2 &2(ci - c 2 )] 
i(ai6ic 2 + aib 2 Ci + a 2 b 1 c 1 - a 2 b 2 c 2 ) 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



(Al) and ( |A3| ) are just trivial extensions of lower order 
inequalities, and (|A5| ) belongs to the set developed by 
Mermin ]^]. The maximal quantum violations, the num- 
ber of the first inequality of each of the 5 orbits, and the 
sizes of the respective orbits are stated in the following 
table: 



meq. 


| orbit | 


qm.viol. 





16 


1 


1 


128 


5/3 


3 


48 


v/2 


6 


48 


V2 


23 


16 


2 



(Al) 
(A2) 
(A3) 
(A4) 
(A5) 



2. Inequalities for four sites 

For n = 4 we just give the number of the first inequal- 
ity of each of the 39 orbits, its size, and the respective 
maximal quantum violations. The index p labels orbits 
including an element which is invariant under permuta- 
tions of the subsy stem s, an d / indicates factorizing Bell 
polynomials (like (|AlD and (A3) for tripartite systems): 



meq. 


| orbit | 


qm.viol. 




OZ 


1 


1 

t P 


^1 9 


1 843 

± . 04:0 


f 


1 094 


Of O 


6 


1536 


5/3 


7 


3072 


1.932 


15f 
xu / 


192 


v/2 

V 


22 


2048 


1.932 


23 


1024 


V5 


24 


1024 


2 


25 


6144 


V3 


27 


3072 


V3 


30 


3072 


v/3 


60/ 


384 


v/2 


105 


128 


v/2 


278 p 


256 


v/5 


279 p 


512 


2.556 


280 


3072 


2.139 


281 


1536 


1.819 


282 


3072 


1.819 



Zoo 


fil A A 

014:4: 


9 0791 


zoo 


1000 


9 078 


987 

ZO ( 


ItJOU 


9 ^9R 


300 


3072 


2 


301 


6144 


5/3 


303 


3072 


1.819 


317 


3072 


2 


318 


1536 


2 


319 


2048 


2.139 


360 


1024 


2.326 


363 


1536 


V3 


367 


1536 


V3 


383p 


256 


2 


831 f 
j 


128 


2 


854/ 


96 


2 


857 


384 


v/2 


874 


384 


2 


1632 


96 


v/2 


1647 


192 


2 


6014 p 


32 


2V2 



The number of the inequality representing the orbit of 
Mermin's inequality is 6014. 
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